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ABSTRACT 

The recent detection of an unusually hard spectral component in GRB941017 
extending to > 200 MeV is hard to explain as a synchrotron emission from 
shock-accelerated electrons. It was argued to imply acceleration of protons to 
ultra-high energy. We show here that the "high energy tail" can be explained 
as emission from shock-accelerated electrons in the early afterglow epoch, taking 
into account the effect of synchrotron self-absorption. High energy observations 
set in this case stringent constraints on model parameters: A lower limit to the 
total explosion energy E ^ 5 x 10 53 erg (assuming spherical symmetry); An upper 
limit to the density of gas surrounding the explosion, n & 10 _2 (S/10 54 erg)cm -3 ; 
A lower limit to the expansion Lorentz factor Tj ^ 200; and An upper limit to 
the fraction of thermal energy carried by the magnetic field behind the shock 
driven into the surrounding medium, 6b j < 10~ 4 . Such constraints can not be 
inferred from keV-MeV data alone. The unusually low value of sbj and the 
unusually high ratio E/n may account for the rareness of GRB941017-type high 
energy tails. Tighter constraints on model parameters may be obtained in the 
future from optical and sub-TeV observations. 

Subject headings: gamma rays:bursts — gamma rays: theory — radiation mech- 
anism: nonthermal 



1. Introduction 

In fireball models of GRBs (Piran 2000; Meszaros 2002; Waxman 2003) the energy 
released by an explosion is converted to kinetic energy of a thin baryonic shell expanding at 
an ultra-relativistic speed. The GRB is most likely produced by internal shocks within the 
expanding shell. At a later stage, the shell impacts on surrounding gas, driving an ultra- 
relativistic shock into the ambient medium. This shock continuously heats fresh gas and 
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accelerates relativistic electrons, which produce synchrotron emission that account for the 
X-ray, optical and radio emission (the "afterglow") following the GRB. The initial interaction 
of fireball ejecta with surrounding gas produces a reverse shock which propagates into and 
decelerates the fireball ejecta (Meszaros, Rees & Papathanassiou 1994). As the reverse 
shock crosses the ejecta, it erases the memory of the initial conditions, and the expansion 
then approaches self-similarity (Blandford & McKee 1976), where the expansion Lorentz 
factor decreases with radius, r BM = (17£ , /167rnm p c 2 ) 1 / 2 r~ 3 / 2 . Here E is the total explosion 
energy (assuming spherical symmetry), n is the number density of the ambient medium. 

The duration T of the stage of transition to self-similar expansion, during which the 
reverse shock "lives", is comparable to the longer of the two time scales set by the initial 
conditions (Waxman 2003): The (observed) GRB duration T GRB and the (observed) time 
T r at which the self-similar Lorentz factor r B M equals the original ejecta Lorentz factor Tj, 
Lbm(7t) = IV Since the characteristic time over which radiation emitted by the fireball at 
radius r is observed by a distant observer is ~ r/4r 2 c (Waxman 1997), T is determined by 
(Waxman 2003) 



where E = 10 53 E 53 erg and n = 0.1n_icm 3 . Note, that the duration is increased by a 
factor 1 + z for a burst at redshift z. 

Observations of GRB941017 show two distinct spectral components (Gonzalez et al. 
2003): A low energy component, with photon energies e 7 ^ 3 MeV, and a high energy com- 
ponent, e 7 ^ 3 MeV. The low energy component shows rapid variability, has a characteristic 
GRB spectrum peaking at e 7 ~ 0.5 MeV, and decays over ~ 100 s. The high energy compo- 
nent has a very hard spectrum, number of photons per unit photon energy dn^/de^ oc e^ 1 , 
and persists over ~ 200 s. The different temporal behavior suggests that the two components 
are produced in different regions of the expanding fireball. The characteristics of the low 
energy component suggest that it is produced by internal shocks, similar to other GRBs. 
The temporal behavior of the high energy component suggests that it is produced during 
the transition to self-similar expansion (Granot & Guetta 2003). 

The hard, <in 7 /cfe 7 oc e" 1 , spectrum is difficult to account for in models where emission 
is dominated by shock accelerated electrons (Gonzalez et al. 2003). This has lead Gonzalez 
et al. to suggest that the high energy tail is due to electromagnetic cascades initiated by 
the interaction of photons with ultra-high energy shock-accelerated protons (Waxman 1995; 
Vietri 1995; Botcher & Dermer 1998; Totani 1998; Waxman & Bahcall 2000). We present 
here an alternative explanation: Electrons accelerated in the forward shock inverse- Compt on 
scatter optical photons, emitted by the reverse shock electrons, to create the observed spectra. 
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A key point, which allows to reproduce the observed hard spectrum, is the modification of 
the synchrotron spectrum by self-absorption in the reverse shock. We show that this effect 
allows to reproduce the observed hard, high energy tail also in the internal shock phase 
(see § 4). However, we consider the latter explanation less likely, due to the weak time 
dependence of the high energy component. 

Granot & Guetta (2003) have recently considered inverse-Compton emission from re- 
verse shock electrons during the transition to self-similarity, as an explanation to the high 
energy tail of GRB941017. They have found that in order for such an explanation to be 
viable, the Lorentz factor associated with fireball expansion should be higher, T > 10 4 , and 
the magnetic field in the fireball plasma should be much lower (well below equipartition), 
compared to values typically inferred from early afterglow observations (see, e.g., Zhang, 
Kobayashi & Meszaros 2003). We show here that the high energy tail may be explained as 
emission from the forward shock electrons, with fireball plasma parameters which are typical 
to those inferred from GRB observations: T ~ 10 2 and magnetic field close to equipartition 
(within the fireball plasma). In this scenario, the inferred density of plasma surrounding the 
fireball is lower than inferred for most other GRBs (the magnetic field strength in the shock 
driven by the fireball into the surrounding gas is poorly constrained by current observations). 
We consider the latter scenario more likely, since it requires a modification of the parameters 
of the environment external to the fireball, rather than modifications of the fireball physics. 
Our analysis further improves on that of Granot & Guetta (2003) in including the effects of 
self-absorption, and in carrying out detailed (numerical) calculations of the spectra, which 
are necessary given the inferred parameter range (see below). 

This paper is organized as follows. In § 2 we briefly discuss the dynamics of transition to 
self-similar expansion and the plasma conditions during the transition. In § 3 we analytically 
derive the constraints that should be satisfied by model parameters in order to allow an 
explanation of the observed high energy tail as inverse-Compton emission from forward shock 
electrons during this phase. In § 4 we present the results of detailed numerical calculations 
of the spectrum, which demonstrate that the observed spectrum may be reproduced when 
the constraints derived in § 3 are satisfied. Such calculations are necessary since the spectral 
shape near the inverse-Compton up-scattered self-absorption frequency is not well described 
by simple power-law approximations. Our conclusions are summarized in § 5. 
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2. The transition to self-similar expansion 

During the transition, plasma shocked by the reverse shock expands with Lorentz factor 
close to that given by the self-similar solution, 

r^i4o (^) 1/8 ^ 3/8 , (2) 

where T = 10 2 T 2 s. The unshocked fireball ejecta propagate with the original expansion 
Lorentz factor, Tj > T. The forward shock which propagates into the ambient medium is 
highly relativistic, propagating with Lorentz factor pa Tbm, while the reverse shock is mildly 
relativistic, propagating into the ejecta with Lorentz factor wl + Vi/VsM- 

The radiation luminosity and spectrum are determined, for given E, n and Tj, by the 
fraction e# (e e ) of shock thermal energy carried by magnetic field (relativistic electrons), 
and by the shape of the electron distribution function, which is commonly assumed to be a 
power law of index p = — dlnn e /dln7 e for electron Lorentz factors 7 e exceeding 7 em . The 
values of the parameters €b and e e can not be determined from first principles, since the 
micro-physics of collisionless shocks which determines them is not fully understood. We 
therefore leave them as free parameters in the following analysis, and derive the constraints 
imposed on them by observations. Since the forward shock is highly relativistic, while the 
reverse shock is only mildly relativistic, we allow different values for es and e e in the two 
shocks. Thus, for example, r y e m,r = e e ,r(^p/^e)(rj/rBM)(p — 2)/(p — 1) for the reverse shock, 
and Tem,/ — £ej(jn P / 'n^e^BMip — 2)/(p — 1) for the forward shock (for p > 2). 

Synchrotron emission from electrons accelerated in the reverse shock typically peaks in 
the optical/IR band (Meszaros & Rees 1997; Sari & Piran 1999), 

<RS « J5 = °' 2 n -? r l-Ar,-ltZ,-2 eV - ( 3 ) 

Here, Tj = 10 2 ' 5 r ii2 . 5 , e &tT = 10 _1 e e)rj _i, and €B, r = lCT 2 eB, r ,-2- This energy is similar to the 
photon energy e ssa at which the optical depth for synchrotron self-absorption in the reverse 
shock equals unity, 

ob. _ n o F l/3 l/3 r 2/3 r l/3 2/3 1/3 v U) 

Denoting by 7 c r the Lorentz factor of reverse shock electrons for which the synchrotron 
cooling time is comparable to the ejecta expansion time, and assuming that electrons with 
7 > 7c,r emit nearly 100% of their energy during the dynamical time, the photon energy 
density is (for p — 3) 

^ n 7em,r n _ i A -2 zttV 2 3/2^-1/2-p 2 -3 /r\ 

u p h,s y n,r~Z e e , r u = 2.5 x 10 E<& n_\ T 2 r ij2 .5e e ,_ieB )r _2ergcm , (5) 

Jc,r 

where the thermal energy density is u — 4r BM nm p c 2 . 
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3. Constraints on model parameters 



Electrons accelerated in the forward shock, with Lorentz factors 7 em ,/ ~ ee,/(^p/^e)rBM/2 ~ 
10 4 ' 5 , inverse- Compt on scatter the optical photons, boosting their energy by ~ l1 m f to 
~ 1 GeV. The spectral shape of inverse- Co mpton scattered photons at observed energy 
range, £ 7 < 1 GeV, is similar to the spectral shape of the reverse shock synchrotron spec- 
trum at e 7 ^ 1 eV, as the energy of all photons is boosted by roughly the same factor 
~ 7em,/- However, the spectral shape is expected to smoothen somewhat, because Compton 
scattering is not exactly monochromatic. The low energy synchrotron spectrum is affected 
by self- absorption (Eq. 4). The optical depth increases as photon energy decreases, leading 
to a thermal spectrum dn 1 /de 1 oc at e 7 <C e ssa . This spectrum is harder than observed. 
However, it becomes softer as e 7 approaches e ssa ~ £ 7 ,rs, which may allow to produce a 
spectrum consistent with observations. 

Thus, in order to explain the observed spectrum as due to inverse- Compton scattering 
of self-absorbed synchrotron photons, the following conditions must be met: 

1. The energy of synchrotron photons for which the reverse shock optical depth is unity, 
e 7 = e ssa , should be increased via inverse-Compton scattering by the lowest energy 
electrons in the forward shock, 7 e = 7 em j, to « 200 MeV, i.e. 



2. The lowest energy electrons in the forward shock, 7 e = 7 em j, should lose energy (via 
synchrotron and inverse-Compton scattering) on a time scale longer then T. If this 
condition is not met, cooling of the electrons would lead to a spectrum d lnn e /d hi7 e ps 
—2 at 7 e < 7e m ,/, which will produce an inverse-Compton spectrum dn 1 /de 1 oc £ 7 3 ^ 2 
below 1 GeV, softer than observed. This condition can be written in the form 



where 7/c,c,/ (7c,/) is the Lorentz factor of forward shocked electrons for which the cool- 
ing time due to inverse-Compton (synchrotron) emission is comparable to the expansion 
time. The inverse-Compton cooling time is estimated via t coo ijc = t coo ^ syn {uBj /u p h, sy n,r) , 
where ubj = £bju is the energy density in the magnetic field at the forward shock. 

3. Synchrotron emission from the forward shock electrons, as well as inverse-Compton 
emission from reverse shock electrons, should not exceed the observed flux at the 
10 KeV to 1 MeV range. These requirements can be met by requiring either (a) The 



7em,/ £ ssa 




(6) 



7c,/ ^ 7em,/5 1lC,c,f ^ 7em,/5 



(7) 
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flux to be lower than observed, or (b) The photons to be emitted below 10 KeV. For 
synchrotron emission from the forward shock this implies either 

LlCJ_ = Uph^ > 1Q3) or £ o6 Fg < 1KeV> (g) 
L syn ,f U B ,f 

For inverse-Compton emission from the reverse shock it implies either 

7^ = T^^llrnjTTJ > ^ , Or £ £ C>RS = S^llm,r < 1 KeV. (9) 
-^IC,r ^e,r lem,r 

A lengthy, but straight forward, manipulation of the above constraints lead to the 
following constraints on model parameters 1 : 



a. 


E53 


> 


4T 2 3/2 , 


b. 


(fit) 


> 


50 T 2 4 , 








c. 


e e,r,-l 


> 




d. 


e B,r,-2 


< 


15T 2 37/8 , 



(10) 

e. e e j ~ 0.33, 

/. e B ,/,-2 < (2.5 x 10- 4 , 5 x 10~ 4 T 2 3/4 , 5 x lO" 3 T 2 1/2 , 1.2 x 10~ 2 T 2 5/4 ) , 

g. r i)2 . 5 > o.76T 2 1/8 . 

It is sufficient to meet 1 of the 4 conditions (f). It should be pointed out here that the 
redshift of the source of GRB941017 is considered a free parameter in the above analysis. 
A particular choice of model parameters implies a particular choice of redshift given the 
measured flux, as illustrated in few examples in figure 1. 



4. Numerical results 

In order to obtain an accurate description of the spectrum, we have carried detailed 
numerical calculations of the emission of radiation from a fireball in transition to self-similar 
behavior. Our numerical model contains full description of cyclo-synchrotron emission, syn- 
chrotron self absorption, inverse- as well as direct Compton scattering, pair production and 
annihilation and cascade processes occurring at high energies. Full description of the model 
appears in Pe'er & Waxman (2004). In the present version, photons emitted after accelera- 
tion of electrons at the reverse shock, serve as seed photons for physical processes occurring 
after acceleration of electrons in the forward shock. 



A detailed description of the derivation may be found at http://www.weizmann.ac.il/^asaf/ 
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Figure 1 presents three examples of the results of our calculations for three different 
sets of parameter values, consistent with the constraints derived in § 3. It demonstrates that 
the process described above may indeed account for the observed hard spectrum of the high 
energy component. 

An alternative explanation of the spectrum, emission from internal collisions within 
the expanding fireball, is also presented. For this scenario, we assumed e# and e e to be 
similar in the two shock waves created as two plasma shells collide. The observed spectrum 
may be reproduced adopting rapid variability, At var 10~ 5 s, and very low value of the 
equipartition fraction, e# ~ 10~ 7 . Here too, the hard high energy tail is due to inverse- 
Compton scattering of self-absorbed synchrotron spectrum. Accounting for ~ 1 MeV GRB 
emission, like the variable low energy component observed at the first ~ 100 s in GRB941017, 
typically requires in the internal shock scenario equipartition fraction close to unity. This, 
and the temporal characteristics of the high energy component, suggest that the high energy 
tail is produced during the transition to self-similar expansion, rather than in internal shocks. 

The following point should be noted here. Gonzalez et al. (2003) present spectra 
averaged over 5 distinct time intervals. The flux carried by the high energy tail shows 
variations between different time intervals, while large uncertainties in the determination of 
the spectral index in individual time intervals make it difficult to determine whether the 
spectral index is time dependent. In principle, time dependence of the high energy flux and 
spectrum provides additional constraints on the model. Such variations would arise in the 
model due to the expansion of the plasma during the transition phase, and due to deviations 
from homogeneity in the fireball shell as well as in the surrounding medium. Thus, the 
observed time dependence may allow to constrain the distributions of density and Lorentz 
factor within the fireball ejecta, and the distribution of ambient medium density. Given the 
relatively large uncertainties in the determination of the time dependence of the flux and of 
the spectral index, we find that such analysis is not warranted by the current data. 



5. Discussion 

We have shown that the high energy spectral component of GRB941017 may be ex- 
plained as inverse- Compton emission from electrons accelerated by the shock driven by the 
fireball into its surrounding medium during the transition to self-similar expansion (Fig- 
ure 1). The high energy tail may also be explained as inverse-Compton emission from 
electrons accelerated in internal shocks (Figure 1), and in the reverse shock during transi- 
tion to self-similar expansion (Granot & Guetta 2003). As explained in the introduction, 
we consider the latter explanations less likely, since they require modifications of the fireball 
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physics (in particular, magnetic field well below equipartition in the fireball plasma), and 
since the temporal behavior is hard to account for in the internal shock scenario. 

We have shown that the high energy data, when interpreted as emission from forward 
shock electrons, allow to put stringent constraints (Eq. 10) on fireball parameters, such as the 
total (isotropic equivalent) energy E, initial Lorenz factor Tj and ambient medium density 
n. Such constraints can not be inferred from keV — MeV observations alone. The inferred 
values of E and Tj are similar to those typical to cosmological GRBs, while the inferred 
constraint E/n > 10 56 ergcm 3 implies a ratio E/n which is higher than typically obtained, 
E/n ~ 10 54 ergcm 3 (Bloom, Frail & Kulkarni 2003). The inferred values of e e j, €B, r and e e ^ 
are consistent with those inferred from other GRB and afterglow observations. The value 
of €bj is usually less well constrained by observations, and is required to be well below 
equipartition in our case. This, and the large ratio of E/n may account for the rareness of 
GRB941017-type high energy tails. 

The unknown redshift of the source of GRB941017 is treated in our analysis as a free 
parameter: A choice of model parameters implies a choice of redshift (see figure 1). Values of 
E close to the lower limit, ~ 10 54 erg, imply low redshift, z ~ 0.1. This is a direct consequence 
of the fact that the gamma-ray fluence of GRB941017 is unusually high, rather than of the 
fact that its high energy spectrum is unusually hard. The fluence, 6.5 x 10 _4 erg/cm 2 , implies 
an (isotropic equivalent) gamma-ray energy release of ~ 10 55 erg for z = 1. 

Figure 1 demonstrates that different scenarios accounting for the high energy component 
of GRB941017, as well as different model parameters in a given scenario, lead to different 
model predictions for the fluxes at optical, X-ray and sub-TeV energy bands. The predicted 
fluxes are well within the detection capabilities of the SWIFT (in optical) and GLAST (at 
10-100GeV) satellites, and of sub-TeV ground based Cerenkov telescopes (e.g. HESS 2 , 
MAGIC 3 , MILAGRO (McEnery 2003), VERITAS (Weekes et al. 2002)). Therefore, 
optical and sub-TeV observations on minute time scale will allow to put more stringent 
constraints on explosion parameters than those given in Eq. 10. 

We thank J. Granot & D. Guetta for helpful discussions. This work was supported in 
part by a Minerva grant and by a grant from the Rosa and Emilio Segre fund. 
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Fig. 1. — Points with error bars are the observed spectrum of GRB941017, averaged over a 
time interval of 100 s to 200 s following the burst trigger (Gonzalez et al. 2003). Detailed 
calculations are presented for three different sets of model parameters, which are consistent 
with the constraints of Eq. 10: {E = 3 x 10 54 erg, n = 0.1 cm" 3 , r< = 300, e B , r = 10 _1 , e B j = 
I0~ 6 ,z = 0.1} (solid), {E = 10 55 erg,n = 0.1 cm -3 , Tj = 200, e B , r = 10~ 3 ,eB,/ = 10" 5 , z = 
0.25} (dashed) and {E = 10 54 erg,n = 0.03 cm -3 , r« = 220,e Bjr = 0.2, e B , f = 3 x 10^ 6 ,^ = 
0.06} (dash-dotted). The dotted line presents an alternative explanation: rapid variability 
in prompt emission, with model parameters: L = 3 x 10 52 erg s _1 , At var = 10 _5 s,rj = 
1500, e B:r = e BJ = 10" 7 , z = 0.15. 



